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Abstract 

We use the mathematical framework of loop quantum gravity (LQG) to study the quantization of 
three dimensional (Riemannian) gravity with positive cosmological constant (A > 0). We show 
that the usual regularization techniques (successful in the A = case and widely applied in four 
dimensional LQG) lead to a deformation of the classical constraint algebra (or anomaly) proportional 
to the local strength of the curvature squared. We argue that this is an unavoidable consequence of 
the non-local nature of generalized connections. 
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I. INTRODUCTION 



Three dimensional quantum gravity can be defined from a number of different points of view. Tlie first of tliese 
was tfie Ponzano-Regge [l[ model of quantum gravity on a triangulated 3-manifold which provides a quantization of 
Regge calculus. The Ponzano-Regge model is a state sum model for 3-dimensional euclidean quantum gravity without 
cosmological constant using the Lie group SU{2). Let M be a triangulated compact 3-manifold, a state of the model 
is an assignment of an irreducible representation of SU (2) to each edge of the triangulation. For each state there is a 
certain weight, a real number. The weight is given by the local formula: 

w= n (-i)'-'(2j+i) n i-iy^+^^+^^ n {u^^]' 

interior edges interior triangles tetrahedra 

where the weight for a tetrahedron is a 6j-symbol and for each triangle the admissibility conditions are given by the 
requirement that ji + j2 + js G Z, while ji + j2 — js > 0, ji + js — j2 > 0, and j2 + ja — ji > 0. The partition function, 
or state sum, is obtained by summing over all possible values of the spin on every edge in the interior of the manifold, 
subject to fixed values on the boundary: 

Z= (2) 

Since the set of irreducible representations of SU(2) is infinite, for some triangulations this gives a finite sum, whereas 
for some other triangulations this is a divergent infinite sum ^ . 

In the late eighties, early nineties mathematicians were looking for manifold invariants with the hope that this could 
help to classify them and this led Turaev and Viro [2] to the definition of a state sum which, under many aspects, 
was very similar to that of Ponzano and Regge. The two main differences were first that Turaev and Viro explicitly 
showed their model to be triangulation independent and secondly they replaced the Lie group SU{2) with its quantum 
deformation UqSL{2). When the deformation parameter q is a root of unity, then there are only a finite number of 
irreducible representations, which means that the edge lengths are not summed up to infinite values, and the partition 
function is always well-defined. A very important consequence of this is that the answer obtained is finite, and so the 
model appears to be a regularized version of the Ponzano-Regge model. 

In particular, for q = e^, with integer r > 3, the weight Wq for a given assignment of spins is given by a formula 
analogous to that for Ponzano-Regge in which every factor depends on q. Crucially, the spins are limited to the range 
0<j<(r — 2)/2, and there is an extra admissibility condition 

ji +j2+j3<r-2. (3) 

The 6j-symbol is replaced by a quantum 6j-symbol, and the factor for each edge by a quantum dimension. The 
Turaev- Viro model is defined by the partition function: 

(r-2)/2 

Z = N-^ E (4) 

where Nq is a constant depending on q, and v is the number of internal vertices. Since the sum is finite, this is always 
well-defined. Moreover, the partition function Zq depends only on the boundary triangulation, the boundary data, 
and the topology of the manifold. 

The natural question is then in what way the Turaev- Viro state sum is connected to Quantum Gravity? The answer, 
first exhibited by Witten Q, and then rigorously proven by other authors was that the Turaev- Viro state 

sum is equivalent to a Feynman path integral with the Chern-Simons action for SU(2)^. (g) SU{2)_^., where k is the 
level, and then the connection with gravity follows from the fact that the Chern-Simons action for this group product 
is related to the Einstein-Hilbert action for gravity with cosmological constant A (Q, 0) if fc^ = An'^/A. Quantum 
groups also enter the quantization of Chern-Simons theory in the so-called combinatorial quantization approach ([9]), 
where a quantum deformation of the structure group is introduced as an intermediate regularization. 

From the perspective of loop quantum qravity (LQG) only the case of vanishing cosmological constant is clearly 
understood. The quantization is in this case a direct implementation of Dirac's quantization program for gauge 



^ One can show that this divergences are due to infinite contributions of pure gauge modes and eliminate them by appropriate gauge 
fixing. This procedure boils down to working with those triangulations where infinite sums are not present [1C|| . 
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systems. The basic unconstrained phase space variables are represented as operators in an auxihary Hilbert space (or 
kinematical Hilbert space Hkin spanned by spin network states) where the constraints are represented by regularized 
quantum operators. A nice feature of the regularization (which is both natural but also unavoidable in the context of 
LQG) is that it leads to regulated quantum constraints satisfying the appropriate quantum constraint algebra. There 
is no anomaly. This feature together with the background independent nature of the whole treatment allows for the 
removal of regulators and the definition of the physical Hilbert space and complete set of gauge invariant observables. 

We give some additional details on the construction of the physical Hilbert space in the A = case, as they 
will allow us to relate the results of LQG in this simple case with what we have mentioned above concerning other 
approaches. The physical inner product and the physical Hilbert space Hphys of 2 + 1 gravity with A = can be 
defined by introducing a reg ularization of the formal expression for the generalized projection operator into the kernel 



of curvature constraint ([I4I, 

" " Tr[NF{A)] ] . (5) 



P = " J]^ 5{F(A{x)Y ^ j D[N]exp {i J [ 



In [l3| it has been shown how, introducing a regularization as an intermediate step for the quantization, this projector 
can be given a precise definition leading to a rigorous expression for the physical inner product of the theory which 
can be represented as a sum over spin foams whose amplitudes coincide with those of the Ponzano-Regge model. 
Moreover, the divergences mentioned above do not appear in the LQG treatment as the problematic triangulations 
(leading to infinities in the Ponzano-Regge state sum) simply do not contribute to physical transition amplitudes. 

The previous discussion provides motivation to investigate the possible quantization, in the context of LQG, of the 
non vanishing cosmological constant case. If such program could be achieved one would expect to be able to understand 
the Turaev-Viro amplitudes as the physical transition amplitudes or physical inner product between kinematical spin 
network states. Unfortunately, at the moment there is little evidence supporting this expectation. Consequently, 
our work concentrates on the very basic starting point of the Dirac program: the study of the quantization of the 
constraints and their associated constraint algebra. Our results reveals a puzzling feature that, we argue, is proper 
to the nature of the kind of regularizations admitted by the LQG mathematical framework, namely the appearance 
of quantization anomalies. As we will show, the nice interplay between symmetry and the regularizing nature of the 
LQG representation of basic kinematical observables, present in the A = case, is lost in the case of non vanishing 
cosmological constant. Thus, the generic appearance of anomalies in the quantum constraint algebra appears as a 
serious difficulty for the completion of the LQG program for this model. 

The paper is organized as follows: In Section [IT] we start with a brief discussion on the UV problem in LQG. The 
analysis of the constraint quantization and the associated algebra is presented in sections IIIHIIV Al 

II. UV DIVERGENCES, THEIR REGULARIZATION, AND AMBIGUITIES IN LQG 

In standard background dependent QFT in order to make sense of products of operator valued distributions one 
has to provide a regularization prescription. Removing the regulator is a subtle task involving the tuning of certain 
terms in the Lagrangian (counter terms) that ensure finite results when the regulator is removed. However, any of 
these regularization procedures is intrinsically ambiguous. Ambiguities associated to the UV regularization allows 
to classify a theory as a renormalizable QFT if there are finitely many ambiguities which can be fixed by a finite 
number of rcnormalization conditions, i.e., one selects the suitable theory by appropriate tuning of the finite number 
of ambiguity parameters in order to match observations. 

Removing UV divergences by a regularization procedure is intimately related to the appearance of ambiguities 
in the quantum theory. This problem is intrinsic to the formalism of QFT and is also present in the context of 
LQG although in a disguised way In LQG one aims at the canonical quantization of gravity in the connection 

formulation. General relativity (Riemannian in 3d and both Riemannian or Lorentzian in 4d) admits an unconstrained 
phase space formulation in terms of an SU{2) connection A^^ and its non Abelian conjugate electric field Ef. At 
the unconstrained level, the phase space is isomorphic to that of an SU{2) Yang-Mills theory. However, there is 
fundamental difference with Yang-Mills theory in that general relativity contains diffeomorphisms as part of the gauge 
group. This additional structure motivates the introduction of a representation of the phase space basic variables 
as operators in an auxiliary Hilbert space where diffeomorphisms are unitarily represented. This important feature 



^ Even though this problem is a feature of quantum systems with local degrees of freedom it manifests itself also in loop quantum 
cosmology due to the peculiar nature of the model whose definition uses the regulating-structures of the full theory [T5|. 
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is crucial for the implementation of diffeomorphism invariance (which in addition fixes uniquely — up to unitarily 
equivalence — the representation of the basic unconstrained (or kinematic) phase space variables as operators in a 
Hilbert space il7i|). However, a prize to be paid with this choice of quantization is the impossibility of representing 
the connection as a quantum operator but only its holonomy (or parallel transport) along one dimensional paths 
(the so called generalized connection). In this way the basic quantum object is not the configuration variable but a non 
local object constructed from it. As a consequence of this, local quantities entering the definition of the constraints 
(such as the curvature tensor F(A)) do not admit a direct quantum analog and must be replaced by expressions 
written in terms of holonomies or Wilson loops along extended paths. This process eliminates the potential UV 
divergences associated to non linearities in the constraints — as it provides a natural point-split-like regularization of 
local quantities — but generically introduces infinite dimensional ambiguities. As explained in more detail below, the 
reason for this is roughly the same that allows for infinitely many possible ways of defining a lattice action for lattice 
Yang-Mills theory when replacing the Yang-Mills connection by its parallel transport in the lattice action. 

The fact that these non linear quantities require regularization is not surprising as it is a standard feature of 
quantum field theories and the difficulty associated to the definition of products of operator valued distributions at 
the same point. However, what is special of LQG is that its kinematical structure obliges one to use a regularization in 
terms of holonomies as these are the fundamental quantum observables in the auxiliary Hilbert space (the generalized 
connections). Even the connection itself (the classical configuration field) must be 'approximated' by a regulated 
version in terms of holonomies along infinitesimal paths. The fact that this process is ambiguous can be illustrated 
by the following simple example. Consider a one parameter family of contractible loops around the point p G S 
for e > such that the limit e — >■ one has Ofe p. If we denote by We the holonomy of the connection A around the 
loop a^, then the following is an infinite dimensional family of regularizations of the local curvature e°^i^*^(A)|p: 

n/2 ^ . 

/[rW,] = '-^^^^ , (6) 

k=l/2 

for arbitrary coefficients Cj and some positive integer n, and where are the standard generators of the su{2) Lie 
algebra, and Xjid) is the character oi g £ SU{2) in the unitary irreducible representation of spin j. The fact that for 
arbitrary cj's and n these are good regularizations of the curvature component e"''' Fli,(A)\p at the point p comes from 
the fact that for smooth configurations of the connection one has 

f[r^W,]=e^'FUA)\p + ffie'), 

independently of the freedom in choosing the parameters of the family. Moreover, among the parameters of this family 
one should also consider the functional arbitrariness entering in defining the above homotopy of loops around p. This 
is a particular example that illustrates the nature of the ambiguities involved when replacing the local connection by 
its parallel transport. 

This (infinite dimensional) regularization ambiguity associated to the need to replace connections by generalized 
connections are present in the LQG quantization of three dimensional gravity for A = 0. In fact there is an infinite 
dimensional set of quantizations of the constraints all of which are anomaly free at the regulated level. The requirement 
of the quantum constraints to reproduce the appropriate gauge symmetry algebra does not reduce the ambiguities 
in this case. However, this by itself does not necessarily mean that the quantum theory is ill defined: a detailed 
analysis shows that when the regulator is removed in the construction of the generalized projection operator P 
(of equation ^) the physical Hilbert space is uniquely determined: all physical quantities are independent of the 
regularization ambiguities. 

A second source of ambiguities, intimately related to the replacement of connections with (group valued) generalized 
connections, is that the conjugate momenta Ef become non commutative [l6| . This introduces additional ordering 
ambiguities in cases where one is confronted with the quantization of non linear functionals of the electric field. These 
ambiguities are not present in the quantization of three dimensional gravity with vanishing cosmological constant as, 
in that case, all the constraints are at most linear in Ef . However, the ordering ambiguity will arise in the A case 
where constraints are quadratic in Ef . We shall see that using the classical constraints algebra as a guiding principle 
the ordering ambiguities can indeed be reduced (see Section [IV A[) . 

Finally, there is the ambiguity associated to the details of the extended structures used for regulating the basic local 
quantities defining the constraints. For example, the shape of loops used in the definition of Wilson loops used to 
regulate the curvature strength, or the choice of co-dimension one surfaces used when the local E"" field are replaced 
by flux operators admitting quantization in the LQG framework. Although these ambiguities are not really relevant 
in the A = case, they have a strong effect on the constraint algebra in the A 7^ case as well as in more general 
settings. 
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III. CANONICAL THREE DIMENSIONAL GRAVITY WITH A / 



We are interested in (Riemannian) three dimensional gravity with cosmological constant in the first order formalism. 
The space-time is a three dimensional oriented smooth manifold and the action is given by 

S[e,io]= I tr[eAF(cj) + ^eAeAe] (7) 

where e is a su(2) Lie algebra valued 1-form, F{uj) is the curvature of the three dimensional connection uj and Tr 
denotes a Killing form on su(2). Assuming the space-time topology to be = SR where E is a Riemann surface 
of arbitrary genus, the phase space is parametrized by the pull back to S of w and e. In local coordinates we can 
express them in terms of the 2-dimensional connection A]^ and the triad field Sj" — e^'^e^rjjk where a = 1, 2 are space 
coordinate indices and i,j = 1, 2, 3 are su(2) indices and e"*" = — e''" with e^^ — 1 (similarly eab — —^ba with eu — 1). 
The Poisson bracket among these variables is given by 

{Alix),E'^iy)}^SlS}S^'\x^y) . (8) 

Due to the underlying SU{2) and diffeomorphism gauge invariance the phase space variables are not independent and 
satisfy the following set of first class constraints. The first one is the analog of the familiar Gauss law of Yang-Mills 
theory, namely 

G, = DaE^ = , (9) 

where Da is the covariant derivative with respect to the connection A. The constraint (|9]) is called the Gauss constraint. 
It encodes the condition that the connection be torsion-less and it generates infinitesimal SU{2) gauge transformation. 
The second constraint reads 

= e'^\F^,{A)+Ae^^,eiet) = 

= e'^'F^.iA) + Ae^dc'^'^E^Et = , (10) 

where in the first line we have written the constraint in terms of the triad field, while in the second line we have 
used the electric field. This second set of first class constraints is associated to the diffeomorphism invariance of 
three dimensional gravity. In order to exhibit the underlying (infinite dimensional) gauge symmetry Lie algebra it is 
convenient to smear the constraints (jlOp and © with arbitrary test fields a and iV, which we assume not depending 
on the phase space variables ^, they read: 

G{a) = / a'G, = / aWaE^ = (11) 



and 



C{N) = / N^C = / N,{F'{A)+Ae'^''EjEk) = . (12) 

The constraints algebra is then 

{C{N),C{M)} = AG{[N,M]) 
{G{a),G{f3)} = G([a,/3]) 

{C(iV),G(a)} = G([iV,a]), (13) 

where [a, bY — ^^.a^h^ is the commutator of su{2). For future use it will be convenient to split the constraint C{N) 
as 

C{N) = F{N) + E{AN), (14) 

where 

F{N) = / N,{F\A)), E{AN) = / Ae'^'' N,E,Ek. (15) 

JE JE 



In [ISlI one obtains a constraint algebra which mimics the constraint algebra of 4d gravity by using phase space dependent smearing 
functions. This choice considerably complicates the quantization of the (so defined) new constraints as their algebraic structure is no 
longer a Lie algebra. 
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A. Canonical Quantization 



The canonical quantization of the kinematics (i.e. the definition of the auxihary Hilbert space where the constraints 
are to be quantized) is well understood. Details can be found in 19] . Following Dirac's quantization procedure one first 
finds a representation of the basic variables in an auxiliary Hilbert space ,^aux ■ The key ingredient is the background 
independent construction of this auxiliary Hilbert space. The main input is to replace functionals of the connection 
(taken as configuration variables) by functionals of holonomies along paths (the so-called generalized connections) 
7 C S: these are the basic excitations in terms of which the Hilbert space is constructed. Given a connection A and 
a path 7, one defines the holonomy by 

h^\A\ = Pexp / A. (16) 
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The conjugate momentum (densitized triad) £'f field is associated to its flux across co-dimension one surfaces. One 
promotes these basic variables to operators acting on an auxiliary Hilbert space where constraints are to be represented 
by quantum operators satisfying the operator quantum analog of 

The auxiliary Hilbert space is defined by the Cauchy completion of the space of cylindrical functionals Cyl, on 
the space of (generalized) connections A'^ . The space Cyl is defined as follows: any element of Cyl^ 4'rj[A], is a 
functional of A labeled by a finite graph T C S and a continuous function / : 5[/(2)^*'^'"^ — > C where A'^(r) is the 
number of links of the graph F. Such a functional is defined as follows 

^yj\A\^ j{h^,\A\--- ,h^,^^,^\A\) (17) 

where /i-yj^] is the holonomy along the link 7^ of the graph F. If one considers a new graph F' such that F C F', 
then any cylindrical function \['rj[A] defined on F can be promoted to a cylindrical function \I'r'j'[^] defined on F' 
in a natural way [2(?|. Given any two cylindrical functions ^ri,/[^] and ^r2.3[^]j their inner product is defined by 
the Ashtekar-Lewandowski measure 



= / n '^''^■^('^'ri'--- 'S«,(ri,))3('»7i'--- '^7iv,(ri2)) (1^) 
i=l 

where dhi corresponds to the invariant S'C/(2)-Haar measure, F12 C S is a graph containing both Fi and and we 
have used the same notation / (resp. g) to denote the extension of the function / (resp. g) on the graph F12. The 
auxiliary Hilbert space M'aux is defined as the Cauchy completion of Cyl under (|18p . 

The (generalized) connection is quantized by promoting the holonomy (|16p to an operator acting by multiplication 
in ,^aux as follows: 

/I^^-i^] = . (19) 

The triad is associated with operators in ^aux defining the flux of electric field across one dimensional lines. Namely, 
for a one dimensional path r\ and a smearing field a : S — )■ su(2) we define 

E{^,a) EE [ a'E^dx''eab. (20) 



The associated quantum operator in J^aux can be defined from its action on holonomies. More precisely one has 

E^) > h,[A] = I f^lf^ ■ ^ *^^gct of 7) (2^) 
' ' 2 \ h^[A]a (for 77 source of 7) ' ^ ' 

where the curve 7 is assumed to have one of its endpoints at 77, o^^ = ±1 is the sign of the orientation of the pair of 
oriented curves in the order (77,7), and where £p = KG is the Planck length in three dimensions (the action vanishes 
if the curves are tangential to each other). In terms of the triad operator we can construct geometric operators 
corresponding to the area of regions in E or the length of curves. The operators ([19]), and ((2T|) are the basic extended 
variables in terms of which we shall regularize and quantize the constraints (|lll) and ([T^ . So far we have not specified 
the space of graphs that we are considering as this will also be part of the regularization procedure studied in the 
following section. 



* A generalized connection is a map from the set of paths 7 C S to SU{2). It corresponds to an extension of the notion of holonomy 
h-y[A] introduced above. 
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IV. DISCRETIZATION PRESCRIPTION 



According to the canonical quantization program, we now need to introduce a quantization of the smeared con- 
straints (jlip . and (|12p preserving the gauge symmetry algebra at the quantum level. As we mentioned in previous 
sections, before quantizing it is necessary to translate the classical variables entering the definition of the constraints 
in terms of holonomies of the connection and fluxes of the electric field. In order to do that we need to define a discrete 
structure on top of which we can construct these extended variables. We do so by introducing an arbitrary finite 
cellular decomposition of S. We denote n the number of plaquettes (2-cells) which from now on will be denoted by 
the index p e C^. We assume the plaquettes to be squares with edges (1-cells denoted e e Ce) of length e in a local 
coordinate system. It will also be necessary to use the dual complex Cs* with its dual plaquettes p* S Cs* and edges 
e* G Cy,* (see Fig. [2]). Both cellular decompositions inherit the orientation from the orientation of E. The cellular 
decomposition defines the regulating structure. We now need to write the classical constraints in terms of extended 
variables in such a way that the naive continuum limit is satisfied. Namely, it is necessary that (for smooth field 
configurations) the regulated classical constraints become the classical constraints in the limit e — > (or equivalently 
n — )• oo). 

Consequently, the phase space variables Ef and A\ are discretized as follows: the local connection field is now 
replaced by the assignment of group elements h{e) = Pexp(— A) g SU{2) to the set of edges e e C^. We discretize 
the triad field Ef by assigning to each dual 1-cell e* the su(2) element Ei{e*) = J^^eabE^{x)dx°', i.e. the flux of 
electric field across the dual edge e*. With this decomposition of S we can write the regularized versions of the 
constraints (fTT|) and (fT^ as: 

Cla] = J2 tr[aP*G^'*] = (22) 

and 

C"[^] = ti'I^^C'T' (23) 

peCs 

where G^* and are explicitly defined below. 

Finally, we must define the set of states to be considered when studying the regularized constraint algebra. The 
allowed states will be a subset Cyl{C^) C Cyl consisting of all cylindrical functions whose underlying graph is 
contained in the one-skeleton of C^. In other words, the allowed graphs must consist of collections of 1-cells e £ C^. 



A. Regularized constraints algebra 

With the prescription introduced in the previous section we are now ready to compute the discrete version of the 
algebra (IT^ . Let us start with the sub-algebra of the Gauss constraint. Due to the simplicity of the action of SU{2) 
gauge transformations on holonomies it is straightforward to obtain expressions for the regularized constraints that 
are anomaly free as far as it concerns the subalgebra of SU{2) gauge transformations. Thus the Gauss constraints are 
more simply quantized by concentrating on their exponentiated versions. Namely, instead of writing the infinitesimal 
generator of SU (2) gauge transformations as a self-adjoint operator acting on J^aux it is simpler to directly construct 
an unitary operator generating finite SU{2) gauge transformations. This follows from the fact that under a gauge 
transformation g :Yi ^ SU{2) the generalized connection transforms according to 

h^[A]^g^'h^[A]g, (24) 

where gs = g{xs),gt = g{xt) G SU{2) are the value of the gauge transformation at the source and target points of 7 
respectively. Finite SU{2) gauge transformations are then represented by the unitary operator Ucig) whose action 
on '^T.f S Cyl is 

C/G(ff)>*r,^({/i^J) ^*r,/({.g-^/i^.ff..}). (25) 

The requirement that the quantization of the constraints satisfy the quantum counterpart of (jl3p translates into the 
following equations for the unitary generators 

UG{9l)UG{g2) = UG{g29l). (26) 
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The previous equations is satisfied by our definition and hence the canonical quantization of the Gauss constraint 
presented here is anomaly free. Given a : E — >■ su{2) one can compute the infinitesimal generator G(a) from the 
previous line as 

G(a) = -ihjUG{exp{ta))\t=o. (27) 

Now the previous definition translates into the graphical action shown below which can be directly obtained by 
concentrating on a single graph node and formally expanding to first order in a equation (I25p where one has replaced 
g ~ 1 + a. This allows us to define the action of TT[aP* G^*] appearing in (1^^ where we have regularized the Gauss 
constraint with a sum over dual plaquettes. In fact, since each dual plaquette correspond to a single node of the cellular 
decomposition, we can define the action of Tr[aP* G^*] on the node (to which the dual plaquette p* corresponds) which 
is target of four holonomies using the following graphical notation: 

TT[aP*GP*] [> hi (g) h2 (g> hs (g) h4, = (28) 




where we have omitted the index for the smearing field in the figure. Notice that from (|2ip and the last line, the 
action of the Gauss constraint on the node can be interpreted as the flux operator i?(C, a) across an infinitesimal 
circle G (a zero area circle) centered at the node which is reminiscent of the well known geometric interpretation of 
the Gauss law in the Abelian case such as electromagnetism. 

The other cases corresponding to all possible orientations of the four holonomies can be obtained from ([?!]) or from 
((27)) in a similar fashion. The action of the Gauss constraint on arbitrary elements of Gyl{Gs) can be obtained from 
the previous equations by the standard rules of differentiation of functions of finitely many copies of SU{2). 

From the action of the Gauss constraint in a given dual plaquette above it is immediate to see how the Gauss 
constraint sub-algebra is preserved also at the quantum level without anomalies, as we expect from equations (j26p 
and ([27)) . Infact, due to the local action of the Gauss constraint it is sufficient to concentrate on the action of the 
commutator [G"-{a), G"(/3)] on a single node. At the regularized level this means that only the dual plaquette around 
the given node is really relevant and so, among all the terms in the commutator due to the sum over dual plaquettes 
in the definition of G", only one gives a non vanishing contribution, namely: 

[G«(a), G"(/3)] [> * = [Tr[aP*GP*],Tr[pP*GP*]] \>hi®h2®hz®hi^ 

= {{[aP*,pP*] hi) g) h2 g) h3 hi + hi g) i[aP*,/3P*] /la) /is /i4 + 

+ hi®h2® {[aP\pP*] hs) ®hi + hi®h2®h3® {[aP\pP*] hi) 

= G"{[a,/3])t>^ . (29) 

Hence, the quantum Gauss constraints reproduce the correct commutator algebra and are thus anomaly free as far 
as the SU{2) sub algebra is concerned. This property is related to the well known, and extremely useful fact in the 
context of lattice gauge theories, that discretization does not break the Yang- Mills gauge symmetry. 

We shift now the attention to the definition of the regulated constraints associated with ()23p . According to (fH)) 
and (jlSp. we can write the regulated constraint as 

C [N] = F'^iN] + E"" [AN] (30) 

which allows us in what follows to consider the regularization of F[N] and E[AN] separately. The main observation 
that motivates the regularization of F[N] is that the integral defining this first term can be approximated by a Riemann 
sum over plaquette contributions according to F[N] = J^tr[NF{A)] — lime_j.o ^^t'^l-^^^^l- The Riemann sum is 
then the regulated quantity to be promoted to a quantum operator. We only need now to approximate the curvature 
tensor by an expression constructed in terms of holonomies of the connection along suitable paths in such a way 
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FIG. 1: Orientation of plaquette holonomies chosen for the regularization of F[N]. 




FIG. 2: On the left: portion of the cellular decomposition Cs (thin lines) and its dual (thick lines). On the right: the 
edges of are shifted toward the corresponding nodes. The flux operators necessary for the definition of the regularization of 
E[AN] are defined in terms of the latter shifted dual edges. 

that the latter expression converges to the local curvature of the connection A (for smooth configurations of A) in 
the limit e — >■ 0. To this end we define the holonomies £ SU{2) associated to each oriented edge e S C^. In 
terms of this definition the holonomy around a single plaquette p becomes = g^p ■ ■ ■ g^p where here 

for i — 1, • • • ,4 are the corresponding edges bounding the plaquette of interest. Finally using that (for smooth 
configurations) = 1 + + 0{e'^) a natural candidate for regulated F[N] is 

F^IN] = ^F[iVP] = ^tr[7Vf(W^f)]. (31) 
p p 

Notice that the previous regularization corresponds to the choice n = 1 at each plaquette in the general formula given 
in equation At the end of this paper we will discuss the generalization of our analysis to arbitrary regularizations. 
We will argue that the main result of this work is indeed generic and thus independent of the choice made here for 
the sake of simplicity. 

The regularization of E[AN] is more subtle as one needs to replace the classical smooth field E°- by extended flux 
operators which, according to (pij) . act as 'grasping' operators. An important observation, that follows from the action 
of the regulated Gauss constraint, is that the grasping operators must act at the endpoints of the edge holonomies 
in order to avoid inconsistency. Therefore, instead of smearing the i?-field along the edges of the usual dual cellular 
decomposition 77 S C* we need to work with the flux of i?-ficld associated with the shifted edges depicted on the right 
of Figured! By abuse of notation we shall keep denoting the shifted edges 77 S C*. We write the regulated quantity 
corresponding to E[AN] as 

E''[AN] = J2mNpl (32) 
p 

where (concentrating on the shadowed plaquette in Figure [2]) 

E[ANp] = Ae,,feiV; {E{m,r^)E{m,r'')) + Ae^.^N; {E {7^2, r^)E {1^3, r'^ j) 

+ Ke^jkK (S(?73,T^)S(ry4,T'=)) + Ae.jfciV; (^(774, r^)£;(??i, r'^)) , (33) 



where 77^ e C* are the four shifted edges shown in the figure that are dual to the shadowed plaquette p € Ce , and the 
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FIG. 3: On the left: a generic non SU{2) invariant state \l/ in our cellular decomposition. On the right: explicit illustration of 
the orientation of the paths used in the regularization of the constraints. 
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FIG. 4: All possible values of the orientations o^-y of the crossing appearing in equation (|2ip . The edges rj regulating the flux 
operators are the dotted lines while the holonomies are along the continue lines 7. 



operators E{ri,a) are defined in (j2ip . Let us start with the quantum version of 

{C[N],C[M]} = AG{{N,M}) . (34) 

In order to compute the discrete analog of it is sufficient to derive the action of the r.h.s and l.h.s of this 
commutation relation around a single vertex, more precisely around a generic non S'?7(2)-invariant state at a vertex of 
the cellular decomposition, since both sides of the commutator (p4|) vanish for N^, belonging to two plaquettes 
p and p' that do not share a common vertex. Therefore the quantum version of eq. has to be computed for a 
state of the form shown in Fig. [3] and the sums over plaquettes inside the expressions of C"[iV] and C^ [Af] give not 
vanishing contributions only for a finite number of plaquettes, precisely the ones around the given vertex. 

Let us call this non SU{2) invariant state VE" and write the discrete action of the l.h.s of ([M)) on 5': 

[C"[7V],C"'[M]]c>* ^ 

= [[C(7Vi),C(Mi)] + [C{N^),CiM^)] + [C(7V1),C(m3)] + [CiN^),CiM^)] + 

+ [C{N^),CiM^)] + [C(iV2),C(Mi)] + [C{N^),CiM^)] + [C{N^),CiM*)] + 

+ [C(7V3),C(m3)] + [C(iV3),C(Mi)] + [C(iV3),C(M2)] + [C{N^),C{M*)] + 

+ [C{N^), CiM*)] + [C{N*), C{M^)] + [C{N*), C{M^)] + [C{N^), C{M^)]]l>^. (35) 

We have now to choose also a convention of signs for the action of each E field on holonomies, we pick the one shown 
in Fig. [31 where the dashed arrows represent the E field and the black ones holonomies. With these conventions 
we are now ready to compute the action psp . Let us start with the terms involving the commutator of regulated 
quantities at the same plaquette. For the plaquette 1 we have: 



[C{N^),C{M^)]\>^ {[[F[N^],E[AM^]] + [E[AN^], F[M^]]) t> * 
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— (same diagrams switching Ni -H- Mi) , 



where we have used the Leibnitz rule when the flux operators act on products of holonomies, and we have symmetrized 
the local action of the graspings at the node. More precisely, when the product of -E-fields appearing in E'^{AN) act 
at a given point there is, on the one hand, a factor ordering ambiguity associated to the non-commutativity of the flux 
operators, and, on the other hand, an ambiguity associated to the place at which the smearing field A'^ is contracted 
with the Wilson line appearing in the regularization of F'^{N). By symmetrization we mean that every time we find 
this ambiguity wc simply sum over all ordering possibilities with equal weights and divide by their number. This is 
the origin of the weight factors in the previous equation. 

It is easy to see that the first three diagrams in the previous equation give an amplitude that is symmetric under 
A^i ^ Ml and so are canceled. This is a consequence of the symmetrized action of graspings acting at the same 
point discussed above. It is important to notice that if wc did not make this choice of (symmetrized) action we 
would immediately get contributions to the commutator that do not vanish when acting on gauge invariant states! 
The avoidance of this obvious anomaly justifies that above choice which is in addition very natural. Thus only the 
contributions of the last four diagrams remain. The result is: 



[C(Afi),C(Mi)][>* = 



A- 



(tr[{r^ r'^jWi] hi e^jhT'^hi /13 O /i4 - tr[{T^ r^jH^^] e^jhT^hi O /12 O /13 O hi) 



- ti[W^] {{[N\M^]hi) O /i2 O /13 O /i4 + /ii O {[N\M^]h'2) O /13 O /14) , 



(36) 



where we have used {t^,t^} = — 5^*-^ in the last line. An analogous calculation shows that: 



[C{N^), C{M^)] > * = - tr[I^^] {hi (3 {[N', M']h2) /13 /i4 + /ii /12 {[N', M^]/i3) hi) 



(37) 



[C{N^), C(M3)] > * = ^ tv[W^] {hi /i2 {[N\ M^jhs) /i4 + /ii /i2 /13 {[N\M^]h4)) , (38) 



[C{N^), C{M^)] > * = - tT[W^] (([AT*, M^]hi) /i2 /i3 /i4 + /ii /i2 /i3 ([^V , M%4)) . (39) 
8 

We have now to compute the commutators among constraints in neighboring plaquettes. Let us start with 
[C{N^),C{M^)]: 



[C{N^),C{M'^)] > * = {[[F[N\E[AM^]\ + [E[AN\ F[M^]\) > * 
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■ s 


•M2) 













= - | {NlMl tr[{r^T"}V7l] /ii /i2 e«fer^>J3 ® + M^A^i tr[{T^ t"}^^] ^^^.^^fe/j^ ^h^^h:,® /14) 

16 16 

An analogous computation shows that the commutators among constraints in plaquettes which just share a vertex 
((1,3) and (2,4)) vanish, while for the other contributions we have: 



A 



A 



C{N^),C{M^)] [> * = —tr[W^]hi®h2®hs, ® [N^,M'^]hi + —iv[W^]hi ® [N\ M'^]h2 (E) hs (g) h 



16 



16 



[C(7V2), C{M^)] t> * tT[W^] hi (g) [N^, M^] h2(i)h^®hi + ^ tr[M^2] 



[C(7V'),C(Ar)] [> * 



A iY\W^] hi ® ^7V^ M^l «) /13 ® /i4 + 4 trfM/'^l hi (g) h2 (g) hs (g [N^, M^l /14 , 
16 16 



[C(iV3), CiM^)] > * = ^ tr[W^3] [Tv^, M^] hi (g, h2 g) hs (g, h^ + ^ tT[W^] hig)h2g) [A^^ M^] /13 /i4 , 



[C(7V'*), C(Af 3)] > ^ = tr[M^3] ^^4^ ^^3] h^^h2(gh3(ghi + ^ tr[W^4] /ii (g) /12 (g [^^^ M^] /13 (g /i4 , 



A_ 
16 



\C{N^), C(A//i)] > * = A tr[T4^i] hi (g h2 (g h^i (g [iV*, /i4 + A trfVK^i ^ r^4^ ^jii h2(gh:i(ghi. 
16 16 

The fact that four different values of the smearing field enter the action of the regularized constraint at a single point 
is an artifact of the discretization since the values of the smearing field iV, being discretized, changes from plaquette 
to plaquete in the four plaquettes surrounding the node of interest. The 'multi-valuedness' of smearing fields created 
by the discretization makes more obscure the meaning of the previous equations and should not affect the final result. 
A more transparent interpretation of the equations is obtained by coarse graining the discrete smearing fields (the 
equalities that are produced by this means will be denoted by the symbol = , where e.g. stands here for coarse graining) . 
Therefore we can evaluate the action ([55)1 by the appropriate coarse-graining procedure = N"^ = = N'^ = N, 
= M"^ = = Af^ = M and see that at the quantum level the algebra of the curvature constraint with itself 
presents an anomaly, namely: 

[C^[A^],C^[M]] o * = 

c=, 4(tr[W^^] + tr[T4^2] ^ ^^j^sj ^ tT[W^]) [{[N, M]hi) (g h2 g) hs (g h^ + 


+ /ii (g {[N, M]h2) g) hs h4 + hi g) h2 {[N, AfJ/ia) (g /i4 + /ii (g ^2 «) /13 «) {[N, Af]/i4)] = 



-[N,M]) C> * , 



(40) 
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where G« ( ^^f^ [iV, M] ) = T,viJ2pev^AWp]/8){T,pevGP[N, M]). Notice that the classical analog reduces— once the 
regulator is removed (e 0) for smooth field configurations and due to the fact that Tr[M^]/2 — )• 1 — to the correct 
algebra: 

{C[N],C[M]}^AG{{N,M}). 

The presence of tr[VF] in the regulated algebra will not disappear in the quantum case: this is a genuine quantization 
anomaly introduced by our regularization method. 

We can now compute the action on the state 5* of the commutator between the scalar and the Gauss constraints 
and verify that the relation 



[C"(iV),G"'(M)] = CiiN^M] 



(41) 



holds. In this case when acting with the Ihs of (j4T|) on a single node, for the curvature constraint there will be four 
relevant terms in the sum over plaquettes entering the definition of C {N) (the terms related to the four plaquettes 
around the given node), each with a different value of the smearing field N^, while for the Gauss constraint there is 
only one relevant term (the one related to the dual plaquette around the given node) and we will denote M^* the 
value of the smearing field associated to this term. Thereby, for the Ihs of (j4ip there are four commutators to compute 
and for each of them there are two terms: one related to the part of C and one to the E'^ part. 

Since the F{N'''), for all the four plaquettes {p = 1,2,3,4) around the given node, depend only on holonomies the 
Leibnitz rule applied to the their commutators with the Gauss constraint in the dual plaquette around the given node 
G[MP*] implies directly 



[F"(iV),G"'(A/)] [> = F"([iV,A/]) t> 

This implies that we can simply concentrate on the commutator of E"{AN) with the Gauss constraint. 
Starting with the plaquette 1, we have, omitting the index for the smearing field MP*: 

[E{AN^),G{M)][>-i' = 



(42) 
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(43) 



Simple algebra yields 



[£;(A7V^), G{M)] O * = ■J^^'ii Mh2 (g) /i3 (g) /14 - ^Mhi g) N^h2 g) /13 (g hi, 



(44) 



In an analogous way one can compute: 



[E{AN^),G{M)] [> * = jhi (g N^h2 (g Mhs (g - ^/ii g) Mh2 g) N^h^ (g /i4 , 



(45) 
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[E{AN^),G{M)] [> * = jhi (g)h2(E) N^hs (g) Mh^ - jhi /12 M/13 N^h^ , (46) 
[E{AN^), G{M)] [> * = I^M/ii (g) /i2 «) /13 «) ^^/i4 - ^^^^1 <8) /i2 «) /i3 • (47) 



We can now sum up all the contribution and we get: 

[E'^{AN),G''{M)] = 



+ ^N^hi (g) M/12 (g /13 (g - ^M/ii g) 7V^/i2 g) /13 g) + 

At a o 

+— /ii g) 7^^/7,2 Mh3 g) ft,4 - — /ii Mft,2 ® N^hz ft.4 + 

A c A , 

+ -/?,! (g ft,2 «) ^ /i3 «) ^^4 - -^hi g) ft,2 g) M/13 g) N-^hi^ + 

+ ^M/ii g) /i2 «) ^3 A^^^4 - -J^*^! /l2 «) /13 g) M/l4 . (48) 



We have now to compute the r.h.s of (HI]). Starting again from the plaquette l,we get: 



E{A[N^,M]) [> * = ^NlMje.^k r'-h^ g) r'h^ ^h^i^h^ 



^N'^hi g) M/12 g) /13 <8) ft.4 - ^M/ii g) 7V^/i2 g) /13 <8) ft.4 ■ (49) 



An analogous computation shows that: 



£:(A[iV^, Af]) t> * = ^/ii g) iV^/i2 M/13 (g /i4 - ^/ii g) Mft2 «) A^^/i3 (g /14 , (50) 

, A , A o 

£:(A[iV^, M]) t> * = —hi ®h2® N-^hs (g) Mh^ - —hi g) /12 <8) M/13 g) N-^hi , (51) 



£:(A[iV'', Af]) > * = ^ AT/ii g) /i2 g) /13 g) iV''/i4 - ^^^'ii ®h2(gh3® Mh^ . (52) 



Thus, summing all the four contributions to the action of the r.h.s of (|4T|) on 5*, we get: 

£;"(A[Ar,A/]) t> * = 

+ ^N^hi g) Af/l2 g) /l3 g) /).4 - ^Af/li g) iV^/l2 (g /l3 (g /l4 + 

A o A o 

+— /(-I g) iV'^/i2 ® Mhj, g) ft,4 - — ft-i ® A/ft,2 ® N^hs ® ft,4 + 

A o A , 

+— /ii ®h2® N-^hs g) Arft,4 - — /ii ®h2® Mh3 g) 7V''ft,4 + 

+ ^A//ii ig) h2 <S) hs <Si N'^hi - ^N'^hi (g /i2 (g ft.3 g) A//i4 . (53) 

From eq.gll) and dSS]) we conclude that [E''{AN),G'^{M)] = E''{K[N,M]) which combined with ^ yields 

[CiN), G^iM)] = C"([A^, M]), (54) 

as expected. 

To summarize, we have seen that the quantum version of the constraints algebra ([T^ of gravity in 2 + 1 dimensions 
with non-vanishing cosmological constant reads: 

[C-(A^),C-(Af)] = AG-(^^[iV, Af]) 
14 



[G«(iV), CiM)] = CilN, M]) 



[CiN), G"{M)] = C"i[N, M]) (55) 

Relations (|55|) show that just the commutators among the scalar constraints present an anomaly due to the presence 
of the factor Hffl. jn the smearing of the Gauss law on the r.h.s. of the first equation. We see that the regularization 
does not break the internal gauge group SU{2); however, it does break the part of the gauge symmetry group related 
to spacetime diffeomorphisms. Notice also that the anomaly is a genuine quantum effect. If we had computed the 
Poisson algebra of regularized constraints instead we would have found basically the same result (where commutators 
are replaced by Poisson brackets). However, in that case the problematic factor disappears in the continuum limit as 
iim = 1 + ^(e4)_ 

It is useful to rewrite the result in terms of a constraint reparametrization that exhibits more clearly the compact 
character of the gauge symmetry group of Riemannian 2+1 gravity with positive cosmological constant. At the 
classical level we know that the algebra (fT3|) generates a local su(2)® su(2) symmetry, as one can see immediately by 
defining F±(7V) = C{N) ± y/KG{N) and computing: 

{F^{N),F^{M)} = ±2^fKF^{[N,M]) {F^ {N) , F^ {M)} = Q . (56) 

In order to see what modification to this local symmetry appears at the quantum level we define the operators 
F«±(iV) = CiN) ± VAC(iV) and then compute the commutators [F«±(iV), i^«±(M)] and [J^«±(A^), F«=F(Af)]. 
The result is 

[i^«±(A^),F«±(Af)] = ±2\/AF^±([iV, Af]) + :|G"((tr[M^] - 2)[iV,M]) (57) 

and 

[F«±(7V),i^«T(M)] = |G-((tr[W^]-2)[iV,A/]). (58) 
V. DISCUSSION 

We have precisely computed the regulated quantum constraint algebra for Riemannian 2+1 gravity in the connection 
formulation of LQG. The nature of the kinematical Hilbert space of LQG imposes the need of a regularization in 
the definition of the quantum constraints. This is so, due to the fact that only the holonomy (and not the local 
connection) and conjugate fluxes (instead of the local triad) can be quantized in the kinematical LQG representation: 
the fundamental operators representing phase space variables are of extended nature. We have studied a simple 
regularization of the constraints that leads to the correct naive continuum limit of both the constraints and their 
Poisson algebra. However, when these regulated constraints are quantized the regulated quantum constraint algebra 
becomes anomalous due to the presence of plaquette loop operators that do not go away in the refinement limit of 
the regulating lattice ^ . 

There is a large freedom in the choice of regularization of the constraints. One source of ambiguity comes from 
the fact that the flux operators corresponding to the local triad do not commute in the quantum theory. As the 
classical constraints of the theory considered here are non linear in the triad field — due to the presence of a non 
vanishing cosmological constant — this non-commutativity introduces factor ordering ambiguities. We have shown 
that the requirement that the quantum constraint algebra be satisfied completely eliminates this source of ambiguity. 
More precisely, the caculations leading to eq. and finally to eq. (|40p show that if we do not symmetrize over 

all possible orderings then [C{N)^C"\M)] O ^E* 7^ when acting on states 5* annihilated by the quantum Gauss 
constraints. 

Another source of ambiguity can be found in the choice of paths along which the fluxes and holonomies used in 
the regularization are defined. For simplicity we have chosen here a square lattice and its dual. However, it 
should be clear from our calculation that the presence of the anomaly found here is independent of this choice. It 
cannot be removed by playing with this freedom. Finally, there is the ambiguity in the choice of representation used 



^ We would like to point to the existence of results along the lines of this work but in 3+1 dimensions by R. Loll [2 j | where the 3d 
diffeomorphism constraint sub-algebra is tested in a lattice regularization of quantum gravity. 
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in the regularization of the curvature constraint. Instead of regularizing the local curvature using the fundamental 
representation as in (|31|) we could have used the more general function 



= e^N.F' + ^{e^), (59) 

W=l+e'^F 

for arbitrary coefficients cj and some positive integer n. Calculations along the lines presented here show that it is 
not possible to correct the anomaly found here by playing with the parameters of this space of functions. 

Let us mention that the anomaly found here seem to be related to the results obtain in reference [22| by a different 
approach. If we write the constraint algebra as in ([57)1 and (155)) the anomaly is parametrized by the local quantity 
(Tr[iy] — 2). If for a moment we think of this factor as a classical quantity, and we evaluate it on a smooth connection 
field configuartion, then this factor is p roportional to the local curvature squared of the connection. This is in direct 
correspondence with the results of . 

The anomaly found here represents an important obstacle for the implementation of the standard Dirac quantization 
in the LQG representation. In particular the group averaging techniques to solve the first class constraints at the 
quantum level, used in the A = case, are not available. It may be possible to contour this obstacle by means 
of some alternative formulations developed for the quantization of systems with constraint algebras which are not 
associated with a structure Lie group (e.g. Thiemann et al. master constraint program [23]). At this point we can 
only speculate with the possibility that the anomaly found here may be at the end related with the deformation of 
the classical symmetry of gravity leading to the quantum group structure underlying the quantization of 2+1 gravity 
with non-vanishing cosmological constant found by other methods. 



f[NW] 
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